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Introduction 

Microfluidics,  defined  here  as  the  study  of  flows  with  overall  dimensions  ranging  from 
a  few  |um  to  a  few  hundred  |um,  has  been  a  major  area  of  researeh  over  the  last  deeade  [1; 
2],  Diffusion,  which  occurs  over  time  seales  proportional  to  the  square  of  the  length  seale, 
beeomes  mueh  faster  for  flows  at  such  small  length  seales,  and  many  of  the  applieations  of 
mierofiuidies  exploit  this  rapid  diffusion.  Microfiuidie  deviees  sueh  as  “Labs  on  a  Chip” 
(LOC)  and  miero-total  analysis  systems  (pTAS)  are  used  for  high-throughput 
biomoleeular  and  eellular  separation,  reaetion  and  deteetion,  among  other  applieations 

A  major  enabling  technology  for  LOC  is  controlling  the  (usually  ineompressible,  low 
Reynolds  number,  and  internal)  flows  flow  of  “biofiuids” — aqueous  solutions  eontaining 
eleetrolytes  and  charged  biomoleeules — through  nano-  and  mieroehannels.  The  surfaee  of 
the  walls  of  sueh  ehannels  will  beeome  eharged  when  exposed  to  biofluids,  and  this 
(usually  negatively)  eharged  surfaee  will  be  sereened  in  turn  by  a  thin  layer  of  eounterions 
(usually  eations)  drawn  from  the  biofiuid.  This  screening  layer  is  known  as  the  eleetrie 
double  layer  (EDL),  and  it  has  a  eharge  distribution  eharaeterized  by  the  wall  zeta- 
potential  Cw  [3].  The  overall  thiekness  of  the  EDL  is  about  five  times  the  Debye  sereening 
length  A-d  (the  1/e  lengthseale  for  the  potential  from  the  solution  of  the  linearized  Poisson- 
Boltzmann  equation),  whieh  varies  from  a  few  Angstroms  for  aqueous  eleetrolyte 
solutions  at  molar  concentrations  of  0(1  mo  EL)  to  a  few  microns  for  very  pure  water  with 
a  resistivity  of  10  MQ-em. 

In  the  eleetroomostic  (EO)  flows  eommonly  used  in  EOC,  a  voltage  gradient,  or 
external  eleetrie  field  of  magnitude  E,  is  applied  parallel  to  the  wall  to  drive  the  eharged 
fluid  in  the  EDL,  and  the  uniform  flow  of  neutral  bulk  fluid  is  driven  by  viscous  effects  by 
the  nonuniform  “boundary  layer-like”  flow  of  fluid  molecules  and  counterions  in  the  EDL. 
Many  LOC  applieations  employ  EO  flow  beeause  this  uniform  flow  in  the  bulk  gives 
higher  flow  rates  and  less  eonveetive  dispersion  than  the  parabolie  veloeity  profile  (with 
the  same  maximum  veloeity)  of  pressure-driven  Poiseuille  flow. 

Research  Accomplishments  and  Army  Relevance 

The  major  goal  of  this  researeh  is  to  understand  eharged  biomoleeules  in  a  biofiuid 
flowing  through  a  miero-  or  nanoehannel  interaet  with  a  eharged  surfaee  in  the  presence  of 
an  applied  eleetrie  field  parallel  to  the  wall.  In  the  third  and  final  year  of  this  work,  the 
experimental  effort  at  Georgia  Tech  led  by  Dr.  Yoda  has  foeused  on  studying  the 
dynamies  of  eolloidal  partieles,  as  a  model  of  large  eharged  biomoleeules,  in  eombined 
eleetroosmotie  and  Poiseuille  flow  through  mieroehannels,  i.e.,  in  the  presenee  of  an 
applied  eleetrie  field  parallel  to  the  wall  and  flow  shear.  The  objeetives  for  the 
experimental  effort  are  to: 

1)  Determine  how  the  nonlinear  eleetrokinetie  repulsive  lift  foree  reported  by  Kazoe  and 

Yoda  [4]  is  affeeted  by  flow  shear; 


2)  Verify  our  hypothesis  that  tracer  particles  of  different  sizes  in  this  nonuniform  flow 
with  an  electric  field  parallel  to  the  wall  will  have  different  average  wall-normal 
distances,  and  hence  sample  different  average  particle  velocities,  in  this  shear  flow. 

The  objectives  of  the  theoretical  and  computational  effort  at  the  Ohio  State  University  led 
by  Dr.  Conlisk  are  to; 

1)  Calculate  the  velocity  of  a  colloidal  particle  as  it  moves  near  a  wall; 

2)  Predict  the  number  density  of  a  population  of  colloidal  particles  near  the  wall  of  a 
channel; 

3)  Determine  the  origin,  nature  and  magnitude  of  the  forces  on  the  types  of  colloidal 
particles  used  in  the  experiments  in  the  presence  of  shear; 

4)  Determine  using  molecular  simulations,  the  effect  that  different  divalent  cations  have 
on  electroosmotic  flow. 

Experimental  Results  and  Research  Accomplishments 

The  experimental  effort  uses  evanescent-wave  particle  tracking  velocimetry  (PTV), 
which  estimates  velocity  fields  from  the  displacements  of  colloidal  fluorescent  polystyrene 
(PS)  particles  of  radius  a  <  0.5  |um  over  two  successive  exposures  separated  by  a  time 
interval  At.  The  particles,  which  are  (nearly)  density-matched  to  the  fluid,  are  illuminated 
by  evanescent  waves  generated  by  the  total  internal  reflection  (TIR)  of  light  at  the 
refractive-index  interface  between  the  fluid,  an  aqueous  monovalent  electrolyte  solution, 
and  the  fused-silica  wall  of  the  microchannel.  The  method  has  been  used  by  various 
researchers,  including  Dr.  Yoda,  to  study  microchaimel  flows  over  the  first  -0.5  qm  next 
to  the  wall  [5-8].  Although  the  method  is  inherently  limited  to  near-wall  flows,  the  spatial 
resolution  of  evanescent-wave  particle  velocimetry  is  significantly  better  than  the  leading 
velocimetry  method  for  microchaimels,  micro-particle  image  velocimetry  (pPIV),  with  its 
typical  minimum  depth  of  correlation  of  2  pm  [9]. 

In  most  PTV  applications,  the  fluid  velocity  is  assumed  to  be  identical  to  the  particle 
velocity.  In  EO  flow,  however,  the 
particle  velocity  is  not  that  of  the  fluid 
because  the  surface  of  the  particles  (like 
the  walls  of  the  channel)  also  become 
charged  in  the  presence  of  an  aqueous 
electrolyte  solution.  This  particle  surface 
charge  is  characterized  by  the  particle 
zeta-potential  which,  much  like  the 
wall  zeta-potential  is  the  electric 
potential  a  counterion  diameter  away 
from  the  surface  of  the  particle.  The 
external  electric  field  that  drives  the  flow 
will  therefore  also  transport  the  particle  tracers.  For  particle  tracers  with  a  (usually 
negative)  surface  charge  (to  minimize  particle  deposition  upon  the  negatively  charged 
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Figure  1  A  negatively  charged  particle  near  a 
negatively  charged  wall.  Note  that  only  the  mobile 
cations  are  shown  near  the  wall. 


wall),  E  will  drive  the  traeers  towards  the  cathode,  a  phenomenon  known  as 
electrophoresis  [3].  A  tracer  particle  will  therefore  be  convected  by  EO  flow  and  driven 
by  electrophoresis  in  opposite  directions  (Fig.  1).  The  actual  velocity  of  the  particle  Wp  will 
therefore  be  the  superposition  of  the  EO  flow  velocity  Ugo  and  the  particle’s  electrophoretic 
velocity  Mgp,  which  for  the  case  of  “thin”  EDEs  (where  ! aU  1  and  a  is  the  particle 
radius)  is: 


p 

where  s  and  p  are  the  permittivity  and  dynamic  viscosity  of  the  fluid,  respectively.  So  the 
particle  velocity  will  only  be  equal  to  the  fluid  velocity  in  EO  flow  for  ^p  =  0. 

These  particles  are  also  subject  to  forces  along  the  cross-stream  direction,  specifically 
normal  to  the  wall.  It  has  been  known  for  almost  30  years  that  a  «  2.5  pm-5  pm  particles 
convected  by  a  simple  shear  flow,  as  is  the  case  for  Poiseuille  flow  driven  by  a  pressure 
gradient  very  close  to  the  wall,  experience  a  repulsive  force,  known  as  “shear-induced 
electrokinetic  lift,”  due  to  electroviscous  effects  [10].  Dr.  Yoda’s  group  has  also  shown 
that  smaller  a  ~  0.2  pm-0.5  pm  particles  convected  by  a  uniform  flow  in  the  presence  of 
an  electric  field  parallel  to  the  wall,  as  is  the  case  for  EO  flow  driven  by  a  voltage  gradient, 
also  experience  a  repulsive  force,  with  a  magnitude  that  appears,  admittedly  over  a  limited 
range  of  parameters,  to  be  proportional  to  E  and  a  ,  in  agreement  with  theoretical 
predictions  of  a“dielectrophoretic-like”  repulsive  force  due  to  nonuniformities  in  the  local 
electric  field  in  gap  between  the  particle  and  wall  [11].  However,  the  discrepancy  between 
theoretical  predictions  and  experimental  estimates  of  the  magnitudes  of  both  the 
dielectrophoretic-like  repulsion  and  shear-induced  electrokinetic  lift  [12]  are  at  least  an 
order  of  magnitude,  suggesting  that  we  lack  a  theoretical  understanding  of  both  types  of 
wall-normal  forces. 

Given  that  microchannel  flows  are  most  commonly  driven  by  voltage  and  pressure 
gradients,  these  results  suggest  that  the  dynamics  of  tracer  particles  within  1  pm  of  the 
wall  are  quite  complicated,  and  these  dynamics  must  be  accounted  for  to  obtain  accurately 
“map”  particle  velocities  to  flow  velocities  in  EOC.  Improving  our  understanding  of  these 
dynamics  is  therefore  critical  in  improving  the  performance  of  EOC  given  that  particle 
velocimetry  techniques  are  the  leading  technique  for  quantifying  transport  in 
microchannels. 

Evanescent-wave  particle  velocimetry  was  therefore  used  by  Ph.D.  student  Necmettin 
Cevheri  to  study  near-wall  particle  dynamics  of  a  =  125  nm  and  245  nm  PS  particles  in 
combined  EO  and  Poiseuille  flow.  The  major  research  findings  of  this  work  are: 

I.  The  changes  in  the  average  wall-normal  positions  of  the  a  =  245  nm  particles  due  to 
dielectrophoretic-like  repulsion  and  (much  weaker)  shear-induced  electrokinetic  lift 
increase  the  average  velocity  sampled  by  the  tracers 

II.  Unexpectedly  strong  repulsion  is  observed  in  combined  EO  and  Poiseuille  flow,  and 
initial  estimates  of  the  magnitude  of  this  repulsive  force  suggest  that  “the  whole  is 


greater  than  the  sum  of  the  parts” — in  other  words,  the  magnitude  of  this  foree  it  is 
greater  than  the  sum  of  the  dieleetrophoretie-like  repulsion  and  shear-indueed 
eleetrokinetie  lift  foree  magnitudes. 

These  results  are  detailed  next,  and  have  been  presented  at  two  national  and 
international  eonferenees: 

1)  Cevheri,  N.  and  Yoda,  M.  “Nonlinear  eleetrokinetie  repulsion  effeets  in  eombined 
eleetroosmotie  and  Poiseuille  flow  through  mieroehannels,”  65“  Annual  Meeting  of 
the  Ameriean  Physieal  Soeiety  Division  of  Fluid  Dynamies,  San  Diego,  CA  (2012) 
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The  results  in  eonferenee  paper  (item  2)  will  be  submitted  to  Langmuir  by  the  end  of  this 
ealendar  year.  In  addition,  a  paper  summarizing  the  results  from  year  two  of  this  grant  on 
how  divalent  eations  affeet  EO  flow  and  wall  zeta-potentials  was  published  in  early  2013: 

3)  Cevheri,  N.  and  Yoda,  M.  (2013)  “Evaneseent-wave  partiele  veloeimetry 
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Copies  of  this  paper,  as  well  as  the  eonferenee  paper  (item  2)  are  available  at 
https://extranet.aro.army.mil/progress-reports. 

Experimental  Studies 

Experimental  Details 

Eleetroosmotie  and  Poiseuille  flows  were  studied  in  wet-etehed  fused-siliea 
mieroehannels  with  trapezoidal  eross-seetions  of  depth  -'30  pm  and  width  -300  pm  driven 
by  eleetrie  fields  E  =  E'i  (where  E  <  33  V/em  and  i  is  along  the  direetion  of  the 
Poiseuille  flow)  and  pressure  gradients  Ap/L  <1.1  Bar/m,  eorresponding  to  a  simple  shear 
flow  with  a  shear  rate  y<  1600  s  '.  Joule  heating  effeets  were  minimal;  in  all  eases,  the 
temperature  of  the  fluid,  measured  by  a  thermoeouple  at  the  exit  in  some  of  the 
experiments,  was  within  1  °C  of  the  ambient  temperature,  whieh  varied  from  19  °C  to 
21  °C. 

The  working  fluid  was  a  degassed  1  mM  aqueous  sodium  tetraborate  (Na2B407) 
seeded  with  fluoreseent  earboxylate-terminated  polystyrene  (PS)  partieles  of  radii  a  =  125 
nm  with  a  manufaeturer-quoted  polydispersity  of  4.5  nm  (Invitrogen  P881 1)  or  a  =  245  ± 
7.5  nm  (Invitrogen  P88I2).  The  bulk  partiele  number  density  of  the  working  fluid  c^  = 
2.7xl0'^  m^^,  eorresponding  to  volume  fraetions  of  2.2x10  and  1.7x10^^  for  the  a  =  125 
nm  and  a  =  245  nm  partieles,  respeetively. 

The  mieroehannel  was  mounted  on  the  stage  of  an  inverted  epi-fiuoreseenee 
mieroseope  (Eeiea  DMIRE2)  and  illuminated  with  evaneseent  waves  generated  at  the 
bottom  surfaee  of  the  mieroehannel  from  the  total  internal  refleetion  (TIR)  of  an  argon-ion 
laser  beam  eoupled  into  the  substrate  with  prism.  The  intensity-based  penetration  depth  of 


the  evanescent  wave  Zp  =  110  nm-120  nm.  Pairs  of  particle  exposures  with  an  exposure 


time  T  =  0.5  ms  with  a  spacing  within  the  pair  At  =  2  ms  and  a  spacing  between  pairs  of 
200  ms  were  generated  using  an  acousto-optic  modulator  (AOM)  “shutter”  the  laser;  a 
sequence  of  500  image  pairs  (=  1000  images)  were  recorded  for  a  total  data  acquisition 
time  of  about  100  s.  The  fluorescent  particles  suspended  in  the  flow  in  the  center  of  the 
microchannel  (to  minimize  any  effect  of  side  walls)  were  imaged  through  a  63  x 
microscope  objective  onto  an  electron  multiplying  charge-coupled  device  (EMCCD) 
camera  to  give  images  of  the  flow  over  a  region  with  dimensions  of  130  pm  x  37  pm. 

After  identifying  and  removing  the  images  of  overlapping  or  aggregated  tracers,  the 
displacements  of  the  remaining  particles  in  the  plane  parallel  to  the  wall  were  determined 
by  matching  particles  between  the  two  images  in  the  pair  to  their  nearest  neighbor.  The 
particle-wall  separation  h,  or  the  distance  between  the  particle  edge  and  the  wall,  was  next 
determined  from  the  particle  image  intensity  ,  the  area-averaged  integral  of  the 


grayscale  values  in  the  image.  Assuming  that  the  particle  image  intensity  has  the  same 
exponential  decay  as  the  illumination 


I^(h)  =  7°  exp 


(2) 


The  ensemble  of  particle-wall  separations  for  about  3xl0"^  particle  images  was  then  used 
to  calculate  the  particle  number  density  profile  c{h)  over  bins  with  a  width  (dimension 
normal  to  the  wall)  of  20  nm. 


The  particle  displacements  were  divided  into  three  layers,  each  containing  a  similar 
number  of  particle  samples,  based  on  h: 

I.  0<  h<  100  nm 

II.  100  nm  <h<  200  nm 

III.  200  nm  <h<  300  nm 

and  the  particle  velocity  components  parallel  to  the  wall  were  then  simply  the  particle 
displacement  divided  by  At. 


Results  and  Discussion 


Figure  2  shows  the  flow  velocity  profiles  U{z)  for  Poiseuille  flow  (with  no  electric 
field)  on  the  left  at  pressure  gradients  ApIL  =  0.43  Bar/m  {circles),  0.74  Bar/m  {triangles), 
and  1.04  Bar/m  {squares)  estimated  from  a  =  125  nm  {open  symbols)  and  a  =  245  nm 
{filled  symbols)  particles  for  z  <  540  nm,  and  for  combined  EO  and  Poiseuille  flow  at  the 
same  pressure  gradient  values  and  E  =  8.8  V/cm.  For  Poiseuille  flow,  the  data  are 
consistently  higher  than  the  velocity  profile  predicted  by  the  analytical  solution  (solid 
lines);  we  suspects  that  the  discrepancy  may  be  due  to  calibration  errors,  specifically  in 
determining  7°  {cf.  Eq.  2)  since  the  slope  of  the  data  are  in  good  agreement  with  the 

predicted  shear  rate  and  that  the  data  for  both  particle  sizes  are  in  good  agreement.  As 
expected,  the  EO  flow  creates  an  “offset”  in  the  flow  velocity  (proportional  to  E,  over  the 


range  of  electric  fields  studied  here,  though  results  are  not  shown  here),  and  the  slope  of 
the  data  are  in  good  agreement  with  that  of  the  Poiseuille  flow  (dashed  lines). 


Particle  center-wall  distance  z  [nm] 


Particle  center-wall  distance  z  [nm] 


Figure  2  Flow  velocity  profiles  for  Poiseuille  flow  {E  =  0)  at  ISpIL  =  0.43  Bar/m  {circles)  and  1.04  Bar/m 
{squares)  [left]  and  combined  Poiseuille  and  EO  flow  [right]  aiE  =  8.8  V/cm  measured  by  a  =  125  nm  {open) 
and  a  =  245  nm  {filled  symbols)  particles.  The  dashed  lines  represent  the  analytical  solution  for  Poiseuille 
flow  at  the  same  values  of  IspIL.  The  error  bars  denote  the  standard  deviations  over  three  independent 
experiments.  Both  graphs  have  the  same  vertical  axis  label. 


The  distribution  of  such  colloidal  particles  near  the  wall  will  be  nonuniform  because 
electrostatic  repulsion  due  to  particle-wall  EDL  interactions  mediated  by  attractive  van  der 
Waals  forces  due  to  permanent  and  induced  dipole  interactions  [13],  as  described  by  the 
classic  Derjaguin-Landau-Verwey-Overbeek  (DLVO)  theory  of  colloid  science  [3].  If  one 
were  to  assume  instead  that  the  particles  are  uniformly  distributed  over  particle -wall 
separations  0  <  /z  <  300  nm,  the  average  velocity  sampled  by  the  particles  in  this  simple 
shear  flow  should  be  the  flow  velocity  at  the  average  z-position  of  the  particle  center  of  z^ 
=  a  +  150  nm.  Table  I  compares  the  average  velocities  obtained  with  the  particles  of 
different  sizes  with  the  expected  value  of  the  velocity  at  z^  for  the  combined  EO  and 

Poiseuille  flow  case  shown  in  Eig.  1  on  the  right.  Although  particles  of  different  a  of 
course  have  different  velocities  in  this  simple  shear  flow  due  to  excluded  volume  effects, 
the  results  for  fsp/L  =  1.04  Bar/m  demonstrate  that  particles  of  different  sizes  can  measure 
different  velocities  in  the  same  flow — here,  the  velocities  measured  by  the  a  =  125  nm 
particles  are  in  reasonable  agreement  with  the  expected  values,  while  those  measured  by 
the  larger  a  =  245  nm  particles  are  not  because  these  particles  experience  much  stronger 
repulsion  forces,  that  are  discussed  next. 

Shear-induced  electrokinetic  lift  and  the  dielectrophoretic-like  repulsion  were 
quantified  in  terms  of  near-wall  particle  distributions,  expressed  in  terms  of  the  particle 


Table  I.  Comparison  of  average  and  expected  velocities  in  combined  flow  at  is  =  8.8  V/cm  measured  by 
particles  of  different  radii. 
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Figure  3  Particle  number  density  profdes  c{h)  for  EO  flow  (A/»/Z.  »  0)  at  is  =  0  (A),  2.4  (O),  4.7  (O),  7.1 
(□)  and  8.8  V/cm  (V)  [left]  and  for  Poiseuille  flow  {E  =  0)  at  tspIL  =  0.43  (A),  0.74  (•)  and  1.04  Bar/m 
(♦).  In  both  cases,  the  bulk  particle  number  density  c„  =  2.7x10*^  m“^.  Both  graphs  have  the  same  vertical 
axis  label. 

number  density  c  as  a  function  of  the  particle  edge-wall  separation  h.  Figure  3  shows 
c{h)  for  the  a  =  245  nm  particles  in  EO  flow  (with  /SpIL  «  0)  at  E  =  0  V/cm  (A),  2.4 
V/cm  (O),  4.7  V/cm  (O),  7.1  V/cm  (□)  and  8.8  V/cm  (V)  on  the  left,  and  in  Poiseuille 
flow  (with  E  =  0)  at  IspIL  =  0.43  Bar/m  (A),  0.74  Bar/m  (•)  and  1.04  Bar/m  (♦).  For  the 
“baseline”  EO  flow  case  of  E  =  0  (which  is  a  very  weak  Poiseuille  flow)  c  first  increases 
with  h,  reaching  a  maximum  value  at  /t  «  100  nm,  then  decreases  slightly  to  a  value 
comparable  to  the  particle  number  density  in  the  bulk  fluid  =  2.7x10  m  .  The 
number  density  then  decrease  at  a  given  /?  as  E  increases  in  the  EO  flow  cases  on  the  left, 
consistent  with  our  previous  observations  of  the  dielectrophoretic-like  repulsion  force. 
Similarly,  c  decreases  slightly  at  a  given  value  of  h  as  ISpIL  increases  in  the  Poiseuille  flow 


cases  on  the  right,  due  to  shear-induced  eleetrokinetic  lift.  No  such  effects  were  observed 
for  the  smaller  a=  \25  nm  particles. 

Figure  4  shows  c(/z)  for  combined  EO  and  Poiseuille  flow  at  lS.plL  =  0.43  Bar/m  on  the 
left  and  ls.plL  =  1.04  Bar/m  on  the  right,  both  also  at  E  =  0  V/cm  (A),  2.4  V/em  (O),  4.7 
V/em  (0),7.1  V/em  (□),  and  8.8  V/cm  (V).  A  comparison  of  these  number  density 
profiles  with  the  case  of  negligible  shear  shown  in  Figure  3  [left]  shows  that  the  decrease 
in  the  number  of  near-wall  partieles  with  E  is  surprisingly  large  for  these  eombined  flow 
cases,  with  c  values  an  order  of  magnitude  less  than  for  E  >1 .1  V/cm.  This  marked 

decrease  also  appears  to  increase  with  EpIL.  Indeed,  almost  no  particles  were  observed  at 
/z  <  300  nm  at  E  >  8.8  V/cm,  which  is  why  no  cases  are  shown  at  larger  E.  This  marked 
deerease  in  the  number  of  near-wall  particles,  and  the  resultant  increase  in  the  average  z- 
position  of  the  tracers,  is  consistent  with  the  average  velocity  sampled  by  these  particles 
being  signifieantly  larger  than  the  expected  value  at  EpIL  =  1.04  Bar/m  in  Table  I. 


Figure  4  Particle  number  density  profiles  c{h)  for  combined  EO  and  Poiseuille  flow  {ISpIL  »  0)  at  ii  =  0 
(A),  2.4  (O),  4.7  (O),  7.1  (□)  and  8.8  V/cm  (V)  and  ApIL  =  0.43  Bar/m  [left]  and  for  the  same  range  of  is 
values  and  Ap/f  ^  1-04  Bar/m. 

This  large  decrease  also  suggests  that  the  force  driving  these  particles  away  from  the 
wall  likely  exceeds  the  sum  of  the  dielectrophoretic-like  repulsive  and  shear-indueed 
eleetrokinetic  lift  forces.  The  magnitude  of  this  force  was  hence  determined  from  the 
slope  of  the  potential  energy  profile  ^(h)  of  the  particle-wall  interaetion,  which  can  be 
estimated  from  these  measured  c{h)  assuming  that  they  follow  a  Boltzmann  distribution. 
To  try  to  isolate  the  effect  of  the  electric  field  (from  that  of  the  pressure  gradient),  the 
potential  energy  for  the  case  of  “only”  Poiseuille  flow  (z.e.,  when  E  =  0)  was  subtracted 


from  that  combined  Poiseuille  and  EO  flow  at  the  same  ls.pl L.  Table  2  compares  the 
magnitude  of  the  resultant  “additional”  repulsive  foree  (averaged  over  80  nm  <  h  < 
280  nm)  whieh  should  be  due  “only”  to  E,  was  compared  with  the  magnitude  of  the 
dielectrophoretic-like  repulsive  foree  estimated  from  the  cases  shown  in  Figure  3  [left] 
of  EOF  (with  Ap/L  «  0)  at  the  same  value  of  E.  The  results  clearly  show  that  the 
additional  repulsive  foree  F^  is  mueh  greater  than  the  corresponding  value  of 

observed  in  the  absenee  of  the  flow  shear  due  to  Poiseuille  flow.  The  uncertainty  in  these 
foree  estimates  is  quite  large,  about  3  fN  based  on  the  results  of  three  independent 
experiments,  but  the  differences  between  F^  and  F^q  are  greater  than  the  uncertainty  for 

E  >  4.7  V/em.  These  results  in  combined  EO  and  Poiseuille  flow  therefore  show  that  there 
is  ereates  a  strong  wall-normal  foree: 

■  that  repels  negatively  charged  dieleetric  partieles  suspended  in  a  condueting  medium 
with  mobile  eleetrolytes  from  the  surface  of  a  negatively  eharged  wall 

■  with  a  magnitude  that  is  significantly  greater  than  the  sum  of  the  dielectrophoretic-like 
repulsive  force  observed  in  EO  flow  and  the  shear-induced  eleetrokinetic  lift  foree 
observed  in  Poiseuille  flow 

suggesting  that  there  is  a  nonlinear  interaetion  between  the  electrie  field  and  flow  shear. 

Table  II  Estimates  of  the  forces  due  to  the  presence  of  an  electric  field  in  “pure”  EOF  and  in  combined  EO 
and  Poiseuille  flow.  The  uncertainty  in  these  force  estimates  is  about  3  fN,  and  the  estimates  for  at  the 
two  lowest  values  of  is  are  less  than  0.1  fN. 
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Theoretical/Computational  Results  and  Research  Accomplishments 

As  a  first  step  of  the  projeet,  we  reviewed  the  fundamental  theory  of  the  motion  of 
a  eolleetion  of  eolloidal  partieles  near  solid  walls  and  to  eompare  theoretieal 
computations  with  the  experiments  conducted  by  Kazoe  and  Yoda  [4]  for 
electroosmotic  flow.  Understanding  the  near-wall  transport  of  suspended  particles  is 
also  relevant  to  a  number  of  microfluidic  applications  where  a  two-  or  three- 
dimensional  array  of  colloidal  particles  is  assembled  on  a  solid  substrate.  We 
predicted  the  velocity  of  a  single  colloidal  particle  using  Henry’s  Law  of 
electrophoresis  and  found  that  the  presence  of  the  wall  does  not  appreciably  affect  the 
convection  velocity.  In  Year  2  we  developed  methods  of  describing  the  motion  of 
single  colloidal  particles  and  subsequently  show  an  approach  to  describe  the  motion 
and  distribution  of  a  collection  of  particles  (a  mass  transfer  approach)  where  the 
particles  are  considered  to  be  the  third  species  in  an  electrolyte  mixture.  The 
parameter  range  of  interest  in  the  experiments  was  identified  and  the  results  for  the 
particle  velocity  were  compared  with  the  results  of  a  single-particle  electrophoresis 
study.  The  particle  distribution  is  obtained  using  the  aforementioned  mass  transfer 
approach  and  the  results  for  both  the  velocity  and  the  particle  distribution  compare 
well  with  the  experimental  results.  Surprisingly,  the  presence  of  the  walls  has  little  or 
no  effect  on  the  particle  velocity.  When  the  electrostatic  lift  force  on  the  particles  was 
identified,  a  review  of  the  principles  and  fundamentals  of  DLVO  and  non-DLVO 
surface  forces  was  performed. 

In  order  to  analyze  the  corresponding  influence  of  DLVO  and  non-DLVO  forces 
on  the  interaction  of  the  channel  wall  and  a  single  colloidal  particle,  in  Year  3  we 
performed  molecular  dynamics  simulations.  The  results  of  molecular  dynamics  will 
be  compared  against  a  continuum  DLVO  approach.  A  future  objective  of  this  study  is 
to  develop  a  unified  theory  of  the  motion  of  a  distribution  of  colloidal  particles  near 
walls.  We  speculate  that  the  reason  the  particle  velocity  is  independent  of  the  distance 
from  the  surface  as  shown  in  Year  2  work  is  because  the  electrostatic  repulsive  forces 
(EDL  repulsion)  keep  the  particles  far  enough  away  from  the  channel  walls  that 
surface  effects  do  not  affect  the  particle  motion. 

The  force  between  molecules  can  be  purely  physical  in  nature  or  they  may 
involve  chemical  association.  The  van  der  Waals  interaction  should  be  distinguished 
from  the  chemical  one  which  brings  about  molecule  formation  therefore  the  origin  of 
van  der  Waals  forces  is  purely  physical.  In  fact,  van  der  Waals  forces  are  of 
electrostatic  nature,  universal  (acting  on  charged  or  neutral  molecules)  and  non-local. 
They  act  between  all  the  molecules  and  macroscopic  bodies.  Van  der  Waal  forces 
between  the  molecules  depend  on  the  polarizability  of  the  molecules,  i.e.  on  the 
extent  of  the  possible  charge  displacements  inside  the  molecule  or  atoms  as  molecules 
are  not  rigid  entities. 

The  electrostatic  force  between  surfaces  can  be  attractive  or  repulsive;  it  depends 
on  the  surface  charge  density  and  surface  potential  of  each  surface,  and  the  physical 


properties  of  the  medium  separating  the  two  surfaces.  When  surfaces  are  immersed  in 
water,  the  high  dielectric  constant  of  the  water  causes  surface  groups  to  dissociate, 
which  results  in  a  charged  interface.  Furthermore  selective  adsorption  of  ions  from  an 
electrolyte  solution  can  charge  a  surface. 

An  example  of  this  is  the  dissociation  of  silanol  groups  at  the  surface  of  a  silica 
substrate.  The  dissociation  of  the  surface  silanol  groups  creates  a  diffuse  layer  of 
oppositely  charged  ions,  whose  concentration  decreases  as  a  double  exponential  with 
distance  away  from  the  solid- liquid  interface.  When  double  layers  from  opposing 
surfaces  of  equal  charge  overlap,  a  repulsive  force  results  between  the  two  surfaces. 

The  major  research  findings  of  the  theoreticahcomputational  work  are: 

1)  The  velocity  of  a  single  colloidal  particle  using  Henry’s  Law  of  electrophoresis 
and  found  that  the  presence  of  the  wall  does  not  appreciably  affect  the  convection 
velocity; 

2)  Using  a  mass  transfer  approach,  the  concentration  of  the  colloidal  particles  near 
the  wall  is  qualitatively  similar  to  the  experimental  results; 

3)  The  largest  near-wall  force  between  a  colloidal  particle  and  a  charged  wall  is 
likely  the  long-range  electrostatic  force. 

4)  Using  molecular  dynamics  methods,  the  effect  of  divalent  ions  is  tpo  increase  the 
bulk  velocity  in  nanochannels. 

In  the  last  12  months,  the  project  has  produced  one  ISI  journal  article  [25]  already 
published  and  two  papers  in  preparation.  One  abstract  for  oral  presentation  at  the 
annual  meeting  of  the  Division  of  Fluid  Dynamics  (DFD)  of  the  American  Physical 
Society  (APS)  2012  and  at  the  51st  AIAA  Aerospace  Sciences  Meeting  2013.  The 
title  of  the  abstract  and  paper  is  "Effect  of  divalent  ions  on  electroosmotic  transport 
in  a  sodium  chloride  aqueous  solution  confined  in  an  amorphous  silica  nanochanncF’ 
by  Harvey  Zambrano  and  A,  T,  Conlisk.  The  AIAA  paper  will  be  submitted  to  an 
archival  journal  soon. 


Aniurphous  silica  surface 


Electric 

field 


Periodic  huundan 


conditions  in 
X  and  V  directions 


Electrolyte  solution 
»  ith  constant  number  of  ions 


Silica  slabs:  34.76nm  x  2.87nni  x  ~3.10  nm 
Wafer  molecules:  20500 

Ions:  170  inns  (NaC:i  ions  +  divalent  ions  6,  IS  and  .31%) 
Computational  box:  34.76  x  2.87  x  30nm 


Figure  5  Schematic  of  the  system.  The  sketch  shows  the  geometric  specifications  of  the  simulated 
system  for  the  MD  study  of  the  effect  of  divalent  ions. 


Simulation  Details 

Since  details  of  the  simulations  in  Year  1  and  Year  2  are  available  in  previous 
reports,  we  describe  here  only  the  details  of  the  divalent  ion  eomputations. 

Molecular  simulations  with  different  degrees  of  complexity,  have  been 
suecessfully  eondueted  to  study  electrokinetic  transport  in  interfaces,  nanochannels, 
nanopores,  and  earbon  nano  tubes  [16].  In  this  study,  we  perform  Moleeular  Dynamies 
simulations.  We  use  the  MD  paekage  FASTTUBE.  This  paekage  has  been  used  to 
study  several  nanofluidie  systems  including  water  in  and  around  earbon  nanotubes, 
eleetroosmotie  flows  on  silica  substrate  and  capillary  flow  in  silica  nanochannels.  In 
this  work,  we  simulate  a  siliea  nanochannel  as  shown  in  Figure  5.  The  ehannel  has 
dimensions  of  34.76  nm  in  X  direetion,  2.53  nm  in  Y  direetion  and  7  nm  in  Z 
direetion,  being  periodie  in  X  and  Y  direetions.  The  computational  box  is  26  nm  long 
in  Z  direction  in  order  to  avoid  interactions  with  periodic  images  in  that  direction.  A 
SPME  algorithm  with  slab  correction  is  implemented  to  eompute  the  long-range 
electrostatie  interaetions. 

In  order  to  equilibrate  the  eleetrolyte  solution,  MD  simulations  are  performed 
during  2.0  ns  using  a  time  step  of  2  fs.  During  the  simulations  silica  is  kept  frozen. 
Calibrating  the  interaetion  potential  to  mateh  the  experimental  value  of  the 
water  contaet  angle  removes  the  observation  of  ice  [17].  The  system  is  eoupled  to  a 
Berendsen  thermostat  at  300  K  during  the  first  0.5  ns  whieh  is  disconneeted  for  the 
rest  of  the  simulation.  Subsequently,  axial  electrie  fields  are  applied  to  the  systems 
and  NEMD  simulations  are  eondueted  with  a  time  step  of  2  fs  until  the  properties 
(density  and  eleetroosmotie  velocity)  in  the  bulk  region  are  not  ehanging,  during  this 
set  of  simulations  the  system  is  conneeted  again  to  a  Berensend  thermostat  with  a  time 
constant  of  0. 1  ps  to  keep  the  temperature  at  300  K.  We  run  similar  simulations  for 
NaCl  aqueous  solutions  with  different  amounts  of  divalent  ions,  eorresponding  to  6, 
18  and  30%  of  the  total  amount  of  ions  disolved  in  the  solution.  In  all  the 
simulations,  we  keep  electroneutrality  in  the  systems  by  fixing  the  ionic  strength  of 
the  eleetrolyte  solution.  Constant  eleetrie  fields  between  0.2  V/nm  and  1.0  V/nm  are 
applied  to  the  system  in  parallel  direetion  to  the  walls.  The  large  eleetrie  fields 
applied  in  this  study  are  needed  to  enhanee  signal-to-noise  ratio  from  simulations. 
Nevertheless,  the  values  of  the  applied  eleetrie  fields  are  in  line  with  previous  MD 
simulations  of  EOF  in  nanoehannels  [18]  and  with  some  experimental  studies. 

An  aqueous  solution  of  sodium  ehloride  consisting  of  20500  water  molecules  and 
170  ions  is  confined  in  a  amorphous  silica  channel  with  surfaee  charge  density  of 
0.227  e  nm“^.  A  corresponding  Debye  length  of  0.65  nm  is  computed. 

More  than  28  eases  with  different  eonfigurations  are  studied  (two  divalent  ionie 
species  and  4  different  eleetrie  fields  applied  to  eaeh  ease).  For  eaeh  case  a  set  of 
simulations  is  run  during  more  than  40  ns  extracting  data  from  the  last  10  ns.  During 


the  production  run,  the  temperature  of  the  water  will  be  permanently  computed  to 
confirm  the  temperature  of  the  system  is  kept  constant. 

Results  and  Discussion 

Mass  Transport  Model  of  the  Particle  Distribution 

In  this  view,  the  colloidal  particles  are  treated  simply  as  an  additional  species 
characterized  by  its  concentration  within  the  electrolyte  buffer.  The  governing 
equations  for  this  process  are  the  unsteady  Navier-Stokes  equations  along  with 
continuity  and  the  Poisson-Nernst-Planck  system  for  the  electro-  static  part  of  the 
problem. 

In  experimental  studies,  particles  are  generally  assumed  to  follow  the  bulk  fluid 
motion;  however,  it  is  seen  here  that  if  the  particles  are  charged  they  will  deviate  from 
the  bulk  fluid  motion  by  the  magnitude  of  the  electromigration  term,  the  particles  will 
also  deviate  from  the  mean  bulk  fluid  motion  if  diffusion  is  important.  However  this 
term  is  usually  small.  This  method  of  modeling  the  motion  of  particles  has  led  to 
acceptable  agreement  with  experiments  [17]. 

Because  of  the  presence  of  the  bounding  walls  in  a  direction  normal  to  the  walls 
of  the  channel  the  net  flux  of  particles  and  ions  as  well  as  the  solvent  will  be  zero. 
Thus  each  of  the  charged  species  will  have  concentrations  that  follow  the  Boltzman 
distribution. 

The  polystyrene  particles  are  assumed  to  have  a  valence  Zp  that  is  much  greater 
that  the  ionic  species,  and  the  potential  is  calculated  numerically  for  a  variety  of 
valences  for  the  parameters  of  Kazoe  and  Yoda  [4].  The  results  for  the  particle 
number  distributions  are  shown  on  Figure  6.  Note  that  the  computational  results  for 
large  valence  show  a  displacement  from  the  wall  due  to  the  large  negative  value  of 
the  particle  valence.  This  of  course  is  rather  artificial  since  for  such  a  large  valence 
steric  effects  may  arise  that  are  not  included  here.  Nevertheless,  the  simple  model  for 
the  number  density  agrees  qualitatively  with  the  results  from  the  experiments. 
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Figure  6  Particle  number  density  in  EOF  as  a  function  of  the  distance  between  the  wall  and  the 
particle  edge  at  applied  electric  fields  of  a)  0  V/cm;  b)  15  V/cm;  c)  22  V/cm  and  c/)  31  V/cm.  The 
different  line  colors  represent  different  particle  valences. 
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Figure  7  Water  velocity  profiles  for  a  solution  of  sodium  and  chloride  confined  in  a  silica  nanochannel. 
Linear  response  can  be  verified  with  these  profiles  as  different  electric  field  are  applied  parallel  fo  the 
channel  walls. 
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Figure  8.  Water  velocity  profiles  for  a  solution  of  magnesium,  sodium  and  chloride  confined  in  a 
silica  nanochannel.  The  concentration  of  calcium  correspond  to  the  18%  of  the  total  of  the  salt  present 
in  the  electrolyte  solution.  Linear  response  can  be  verified  with  these  profiles  as  different  electric  field 
are  applied  parallel  to  the  channel  walls. 


CoscJ:  27N11  +  91  CI  +  52  Mr 

Figure  9  Water  velocity  profiles  for  a  solution  of  magnesium,  sodium  and  chloride  confined  in  a  silica 
nanochannel  as  a  electric  field  of  0.7  V/nm  is  applied  parallel  to  the  channel  walls.  The  concentration 
of  magnesium  varies  for  each  case  from  0%  to  31%  of  the  total  of  the  ions  present  in  the  electrolyte 
solution.  Higher  concentrations  of  magnesium  result  in  an  electroosmotic  velocity  increase  at  the 
Electric  Double  Layer  (EDL). 

We  compute  the  density  and  velocity  profiles  by  using  the  binning  method  with  a 
binning  size  of  0.12  nm.  Simulations  were  conducted  during  55  ns  and  results  have 
been  extracted  from  the  last  35  ns  of  the  simulation.  Linear  response  of  the  system 
can  be  inferred  from  Figures  8  and  9  that  depict  the  velocity  profile  averaged  over  the 
simulation  length.  Comparing  Figure  7  (no  Mg^^)  with  Figure  5  (with  Mg"^"^),  the 
velocity  is  considerably  higher  in  contrast  to  the  experimental  results  as  described 
above  in  Figure  1.  Moreover,  Figures  8  and  9  show  that  the  velocity  increases  as  the 
concentration  of  the  divalent  ion  increases.  Figure  9  shows  that  the  counter-ions  (Na"^ 
and  Mg^^)  are  attracted  near  the  negatively  charged  silica  walls  while  the  co-ions  CL 
are  repelled  to  the  bulk  fiuid.  However,  Mg^^  ions  are  farther  away  from  the  walls  in 
comparison  to  Na^  ions,  likely  due  to  the  stronger  Mg^^  affinity  for  water.  In  fact,  we 


compute  binding  energies  of  -425  kJ  mol”^  for  the  interaetion  between  Na^  ions  and 
water,  and  of  -510  kJ  mol”'  for  the  interaetion  between  Mg^^  ions  and  water.  The 
eomputed  binding  energies  eonfirm  that  Mg^^  ions  have  a  higher  affinity  for  water 
than  Na^  ions  whieh  results  in  a  more  robust  hydration  in  line  with  the  results  found 
by  Calero  et  al.  [18]. 

These  results  show  that  Mg^^  is  mueh  more  effeetive  dragging  water  moleeules 
due  to  its  hydration  layer  being  mueh  more  robust  so  the  eleetroosmotie  veloeity 
inereases  where  the  proportion  of  Mg^^  is  higher  In  the  bulk  the  loeal 
eleetroneutrality  prevents  ehanges  in  the  eleetroosmotie  veloeity  in  this  area  due  to 
the  proportion  of  eo-ions  and  eounter-ions  keep  the  same  for  different  eases  and 
different  applied  axial  eleetrie  fields.  The  results  suggest  that  by  varying  the 
proportion  of  divalent  ions  in  a  multivalent  eleetrolyte  solution  under 
nanoeonfinement,  it  is  possible  to  exert  eontrol  on  the  eleetroosmotie  veloeity  near 
the  ehannel  walls. 

Summary 

These  surprising  results  suggest  that  the  effeet  of  eombined  EO  and  Poiseuille 
flow  on  eolloidal  particle  dynamies  are  quite  eomplex,  even  though  the  fluid  veloeity 
fields  for  these  ereeping  flows  superpose.  Although  this  project  is  over,  Dr.  Yoda’s 
group  is  still  studying  partiele  dynamies  in  eombined  EO  and  Poiseuille  flow,  at 
present  mainly  for  E  <  0  {i.e.,  when  the  EO  and  Poiseuille  flows  are  in  opposite 
direetions  supported  by  a  grant  from  the  National  Seienee  Eoundation.  Dr.  Conlisk 
will  eontinue  work  on  modeling  and  simulation  of  these  transport  proeesses. 
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